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ABSTRACT 

Multiple carrier-frequency offsets (CFO) arise in a distributed 
antenna system, where data are transmitted simultaneously from 
multiple antennas. In such systems the received signal contains mul- 
tiple CFOs due to mismatch between the local oscillators of trans- 
mitters and receiver. This results in a time-varying rotation of the 
data constellation, which needs to be compensated for at the receiver 
before symbol recovery. This paper proposes a new approach for 
blind CFO estimation and symbol recovery. The received base-band 
signal is over-sampled, and its polyphase components are used to 
formulate a virtual Multiple-Input Multiple- Output (MIMO) prob- 
lem. By applying blind MIMO system estimation techniques, the 
system response is estimated and used to subsequently transform 
the multiple CFOs estimation problem into many independent sin- 
gle CFO estimation problems. Furthermore, an initial estimate of 
the CFO is obtained from the phase of the MIMO system response. 
The Cramer-Rao Lower bound is also derived, and the large sample 
performance of the proposed estimator is compared to the bound. 

feywordi-Multi-user Systems, Distributed Antenna Systems, 
Carrier Frequency Offset, Blind MIMO System Identification 

1. INTRODUCTION 

In any communication system, the received signal is corrupted by 
carrier-frequency offsets (CFOs) due to the Doppler shift and/or lo- 
cal oscillators drift. The CFO causes a frequency shift and a time- 
varying rotation of the data symbols, which need to be compensated 
for at the receiver before symbol recovery. This can be achieved via 
pilot symbols. However, in the case of mobile systems and rich scat- 
tering environments the effects of CFO become time varying and 
even small errors in the CFO estimate tend to cause large data re- 
covery errors. This necessitates transmission of pilots rather often, a 
process that lowers data throughput. In this paper we deal with CFO 
estimation without the need for pilot symbols. In single user sys- 
tems, or in multi-antenna systems in which the transmitters are phys- 
ically connected to the same oscillator, there is only one CFO that 
needs to be estimated. This is typically done via a decision feedback 
Phase Lock Loop (PLL) at the receiver end. The PLL is a closed- 
loop feedback control system that uses knowledge of the transmitted 
constellation to adaptively track both the frequency and phase off- 
set between the equalized signal and the known signal constellation. 
However, depending on the constellation used during transmission, 
the PLL can have an M-fold symmetric ambiguity, resulting in a lim- 
ited CFO acquisition range, i.e., |Ffc| < l/(8Ts) for 4Q AM signals, 
where Ts is the symbol period. Moreover, the PLL typically requires 
a long convergence time. Alternatively, several methods have been 

This work has been supported by NSF under grant Nos. ANI-03-38807, 
CNS-06-25637 and CNS-04-35052, and by the Office of Naval Research un- 
der Grant ONR-N-00014-07-1-0500. 



proposed |2|, |4|, (5), Q (9l that blindly estimate the CFOs and 
recover the transmitted symbols using second-order cyclic statistics 
of the over-sampled received signal. Blind CFO estimation has also 
been studied in the context of orthogonal frequency-division multi- 
plexing (OFDM) systems, where the CFO destroys the orthogonality 
between the carriers (see |3| and the references therein). 

In a spatially distributed antenna system where data are transmit- 
ted simultaneously from multiple antennas, the received signal con- 
tains multiple CFOs, one for each transmit antenna. A PLL does not 
work in that case as there is no single frequency to lock to. The liter- 
ature on estimation of multiple CFOs is rather sparse. Existing liter- 
ature on this topic focuses on pilot based CFO estimation. In | 6|, the 
multiple CFOs were estimated by using pilots that were uncorrelated 
between the different users. In 1 13 1, multiple CFOs were estimated 
via Maximum Likelihood based on specially designed pilots. To ac- 
count for multiple offsets, 1 8 1 proposed that multiple nodes transmit 
the same copy of the data with an artificial delay at each node. The 
resulting system was modeled as a convolutive single-input/single- 
output (SISO) system with time-varying system response caused by 
the multiple CFOs. A minimum mean-square error (MMSE) deci- 
sion feedback equalizer was used to track and equalize the channel 
and to recover the input data. Training symbols were required in 
order to obtain a channel estimate, which was used to initialize the 
equalizer. 

Here we propose an approach for blind identification of multi- 
ple CFOs and subsequent symbol recovery. The received base-band 
signal is over-sampled, and its polyphase components are viewed 
as the outputs of a virtual MIMO system. The time-varying con- 
tribution of the CFOs, together with the transmitted symbols form 
the multiple-input/multiple-output (MIMO) inputs, while the time- 
invariant contribution of the CFOs along with fading channels com- 
prise the system response. By applying blind MIMO system estima- 
tion techniques, the system response is estimated and used to subse- 
quently transform the multiple CFOs estimation problem into many 
independent single CFO estimation problems. Furthermore, an ini- 
tial estimate of the CFO obtained from the phase of the MIMO sys- 
tem response can be used to initialize the PLL eliminating the sym- 
metrical ambiguity problem. The resulting method has full acquisi- 
tion range for normalized CFOs, i.e., \Fk j < P/ {2Ts), where P is 
the over-sampling factor. To evaluate large sample performance we 
establish the Cramer-Rao Bound (CRB) and compare the obtained 
mean square error to it. 

2. SYSTEM MODEL 

We consider a distributed antenna system, where K users trans- 
mit simultaneously to a base station. Narrow-band transmission is 
assumed here, where the channel between any user and the base 
station is frequency non-selective. In addition, quasi-static fading 
is assumed, i.e., the channel gains remain fixed during the packet 



length. The continuous-time base-band received signal y{t) can be 
expressed as 



y{t) = ^ akXk{t - Tk)e 



J 27r Ffct 



(1) 



where represents the effect of channel fading between the fc— th 
user and the base station and also contains the corresponding phase 
offset; Tk is the delay associated with the path between the fc— th user 
and the base station; Fk is the frequency offset of the fc— th user and 
w(t) represents noise; XkXt) denotes the transmitted signal of user 
fc: Xk{t) = Sk{i)p{t — iTs) where Sfe(i) is the i— th symbol of 
user fc; Ts is the symbol period; and p{t) is a pulse function with 
support [0, Tj,]. 

Our objective is to obtain an estimate of s(i) — 



[si(i), 



of the form 



s(i) = AP^s(i) 



(2) 



where P is a column permutation matrix and A is a constant diag- 
onal matrix. These are considered to be trivial ambiguities, and are 
typical in any blind inference problem. 

The received signal y{t) is sampled at rate 1/T = P/Ts, where 
the over-sampling factor P > if is an integer. In order to guar- 
antee that all the users' pulses overlap at the sampling times, the 
over-sampling period should satisfy: Ts/P > rk,k — 1,...K, 
which means that the over-sampling factor P is upper bounded by 
Ts/min{Ti, ...,tk}. Let t = iTs + mT, ra = 1, ...,P ~ 1 de- 
note the sampling times. The over-sampled signal can be expressed 
as 



ym{i) = viiTs + mT) 



- E 



cike 



j27r/fc{iP + m) 



Xk{{i + ^)Ts - Tfc) + wiii + ^)Ts) 



K 



1(3) 



where fk = FkTa/P, (\fk\ < 0.5) is the normalized frequency 
offset between the fc— th user and the base station, and the element 
of the virtual MIMO channel matrix A is given as 



PisTs-Tk), m = l,. 



(4) 



Defining y(i) = [yi{i), ...,yp{i)f ; A = {a„,fc}, a tall matrix 
of dimension P x K; s(i) = [siii)e'^''-'^^'^ , SK{i)e''^''^'''^f ; 



and w(z) = [w{i + j;), w{i + 



eq. ((3) can be written in 



matrix form as 



y{i) = As(i) + w(i) 



(5) 



3. BLIND CHANNEL ESTIMATION AND 
COMPENSATION OF THE CEOS 

Let us make the following assumptions. 

• Al) For each m — 1...P, Wmi-) is a zero-mean Gaussian 
stationary random processes with variance aj,, and is inde- 
pendent of the inputs. 

• A2) For each fc, Sfc(.) are a zero mean, independent identi- 
cally distributed (i.i.d.) stationary with nonzero kurtosis, i.e., 
7*^ = Cum[sk{i) , sl{i) , Sk{i) , sl{i)] / 0. The s^'s are mu- 
tually independent, we can further assume that every user has 
unit transmission power, then Cs = I. 



• A3) The over-sampling factor P is no less than K. 

Under assumption (A2), it is easy to verify that the rotated input 
signals Sk{-) are also zero mean, i.i.d, wide sense stationary with 
nonzero kurtosis. Also, the Sfc(i)'s are mutually independent for 
different fc's. Assumption (A3) guarantees that the virtual MIMO 
channel matrix A in (|5} has full rank with probability one. If the de- 
lays of users are randomly distributed in the interval [0, Ts/P), then 
each row of the channel matrix can be viewed as drawn randomly 
from a continuous distribution, thus the channel matrix has full rank 
with probability one. 

One can apply any blind source separation algorithm (e.g., 1 1 1) 
to obtain 

A = APA (6) 

Subsequently, using a least-squares equalizer we can get an esti- 
mate of the de-coupled signals s(i), within permutation and diagonal 
scalar ambiguities as 



= (A^A)-^A''y(i) = e 



JArg{-A} i 



"'P^s(i) 



(7) 



Denoting by 9k the fc— th diagonal element of Arg{A}, the fc— th 
separated input signal can be expressed as 



Sk{i) = Sk{i)e 



j(-efc+2,r/fciP) 



(8) 



Based on ((Sjl, any single CFO blind estimation method could be ap- 
plied to recover the input signal. Those methods can benefit by a 
CFO estimate provided by the channel matrix estimate as follows. 
The phase of the estimated channel matrix A equals 



* = ArgA = 



V 2nfiP + (t)i 



27V fK + < 



2TTfKP + <t}K 



P (9) 



where (jifc = Ar(;{afc} + Sfe, which accounts for both the phase of 
and the estimated phase ambiguity in The least squares estimate 
of fk is given by 

f 1 nEf=iP*P,fe)-(Ef.iP)(Ep=i*P,fe) 

/fc = t; „,„p — (iU) 



2n 



We can write fk = fk+^k where represents the estimation error. 

Noting that the de-coupled signals Sj{i) in ([8} are shuffled in 
the same manner as the estimated CFOs in dlOb . we can use the es- 
timated CFOs to compensate for the effect of CFO in the decoupled 
signals ([Sj and get estimates of the input signals as 



s(i) 



(11) 



Due to the residual error in the estimated CFOs, we can only 
compensate for a majority of the effect of CFO in ^ and obtain 



Sk{i) = Sk{i)e^^ 



(12) 



Let us apply the PLL to the recovered signals Sj{i) in il2\ . to 
further mitigate the effect of residuary CFO . For 4QAM signals, 
as long as \Pek\ < 1/8, it can be effectively removed by the PLL. 
Thus, the CFO estimator JIOI ) can prevent the symmetric ambiguity 
of the PLL, and can also greatly reduce the convergence time of PLL. 
From ([9}, we can see that the CFO estimator will achieve full acqui- 
sition range for the normalized CFO, i.e., |/fc| < 1/2, which means 
we can deal with all continuous CFOs in the range Fk < P/(2Ts). 



4. CRAMER-RAO LOWER BOUND 

To evaluate the large sample performance of the proposed method, 
we establish the Cramer- Rao lower bound according to f 121. Via 
central limit theory arguments, the received signals y can be approx- 
imated as complex Gaussian signal, with zero mean, and covariance 
matrix given by 



ACsA^ 



AA^ + all 



(13) 



The covariance matrix is valid under assumption Al) and A2). The 
Gaussian assumption of the received signal is reasonable since the 
received signal is a linear mixture of i.i.d. signals. 
Let 



a : 



(14) 



where f — is the vector of unknown CFOs, and 

= [ti, tk]^ is the vector of random delays. The parameter 
we are interested in is the CFOs f , while p and are the nuisance 
parameters. 

Under the previous assumptions and the Gaussian approxima- 
tion, the Fisher Information Matrix (FIM) for the parameter vector 
a is given by 1 12] 



FIM,„=rTr(^C,-^|Sic-^) 

dai " dan 



l,n^ 1,...,2K + 1 



(15) 

Since we are only interested in the CFO parameter f , following the 
derivation in |T2|, we can obtain that 

^CRB"^(f) = G^G - G^A(A^A)"^A"G = G^n^G 

(16) 

where G and A are defined as 



FIM = (^)-(C-®C-)(^) = 



dCy 



A" 



G A 



where Cy = vec(Cy) is a P"^ x 1 vector constructed from columns 
of Cy, and G is of dimension x K, while A is of dimension 
p2 X (K + l). 

To proceed, we just need to evaluate the derivatives of Cy with 
respect to a. First consider dcy/dt^ , it holds that 

TTT = = ^'=^([0 ■ ■ ■ • • • 01A^+A[0 ■ ■ ■ df . . . 0]^) 

dfk dfk 



with dfc 



(afc [1, P]^), where is the Hadamard ma- 



trix product. 

Similarly, we can get dcy/dp^: 



dcy 
dTk 



,dCy, 



with Gfe 



vec([0 ■ 

-p--'^k) 



■0]A"+A[0- 



H 



■Of 



Finally, we have that dcy/dal — vec{Cy ^). Now we have all 
the ingredients to evaluate CRB{f) from l ll6t . 



delays, Tk, k — 1, K are uniformly distributed in the range of 
[0, Ts/P). The input signals used here are 4QAM signals. The es- 
timation results are averaged over 300 independent channels, and 20 
Monte-Carlo runs for each channel. 

The blind source separation algorithm used is 
the JADE method, which was downloaded from 
, http://w ww.tsi. enst.fr~cardoso/guidesepsou.html 

We show the performance of both the pilots-based method and 
the proposed method at different data lengths and SNR set to 30dB. 
For the pilots method, each user transmitted a pilot signal of length 
32, and the pilots were random sequences uncorrected between dif- 
ferent users. In Fig. [T]we show the Mean Squares Error (MSE) 
for the CFO estimator l llOt for different values of the over-sampling 
factor P. To make the comparison fair for different over-sampling 
factor P, the MSE is calculated based on i J^f^^ [{fk - fk)Pf = 

X]fc=i[(^fe ~ Fk)Ts]^- We can see that by increasing P we can 
get more accurate estimates of the CFOs. In Fig. |2l we show the 
Bit Error Rate (BER) for different P's. For both the blind and the 
training based methods, the BER is calculated based on the recov- 
ered signals after the PLL. As expected, the BER performance also 
improves by increasing P. The proposed method appears to work 
well even for short data length. 

Next we show the performance of both methods at various noise 
levels. We set the packet length A*' to 1024. In Fig. [3] we show the 
MSE of the blind CFO estimator dlOl l as well as the training based 
method. We can see that by increasing P we can get more accurate 
estimates of the CFOs. In Fig. |4l we show the BER performance 
after PLL for both blind and training based methods. We can see 
that the proposed blind method has almost the same performance 
to the training based method for SNR lower than 20dB, while the 
training based method can achieve better BER performance for high 
SNR. 

The mean square error of the CFO estimator of JlOt is plotted 
against the stochastic CRB derived in Section|4] In fig. [5] we plotted 
the MSE of the CFOs, as well as the CRB, as a function of the packet 
length T. We can see that the MSE curves are parallel to the CRB, 
and no error floor presented in the plot. Hence there is no bias in the 
estimates and the gap is only due to excess variance in the estimates. 
One possible reason for the existence of the excess variance is that 
we assume that we know the exact channel structure in the derivation 
of the CRB, i.e., the waveform used in transmission, which reduces 
the number of the unknown parameters. In the simulations, however, 
we did not assume any extra knowledge of the channel structure. 

We should note that the PLL is important for good symbol re- 
covery. For example, without the PLL, even if the residual error 
Ptfe = Pfe — Pfe is only 0.001, the constellation will be rotated to a 
wrong position after 0.25/0.001 = 250 samples for 4QAM signals. 
To make sure that the PLL does not have the symmetrical ambiguity, 
we need to guarantee that I Pefc I = \{Fk-Fk)Ts\ < l/8for4QAM 
transmission. Thus, on the average, the maximum tolerable MSE 
for the CFO is in the order of lO^'^. From the simulations, we can 
see that the achieved compensation is sufficient for practical systems 
and commonly used modulation schemes. 



5. SIMULATION RESULTS 



6. CONCLUSION 



In this section, we verify the validity of the proposed method via 
simulations, under the following assumptions. The channel coef- 
ficients flft , k — 1, K are zero-mean Gaussian random vari- 
ables. The waveform p{.) used here is hamming window. The con- 
tinuous CFOs are randomly picked in the range [— 2f^' 2f^)' 



In this paper we have proposed a novel blind approach for identifica- 
tion of a distributed multiuser antenna system with multiple CFOs. 
By over-sampling the received base-band signal, we have converted 
the mulplte-input/single-output (MISO) problem into a MIMO one. 
Blind MIMO system estimation yields the system response, and 



MIMO input recovery yields the decoupled transmitted signals, each 
one containing a CFO. By exploring the structure of the MIMO sys- 
tems response we can obtain a coarse estimate of the CFOs, which 
can be combined with a decision feedback PLL to compensate for 
the CFOs in the decoupled transmitted signals. The proposed blind 
method has full acquisition range for normalized CFOs. We have 
provided a Cramer-Rao bound (CRB) for the proposed blind CFO 
estimators. The analytical results have been validated via simula- 
tions. 
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ABSTRACT 

Multiple carrier-frequency offsets (CFO) arise in a distributed 
antenna system, where data are transmitted simultaneously from 
multiple antennas. In such systems the received signal contains mul- 
tiple CFOs due to mismatch between the local oscillators of trans- 
mitters and receiver. This results in a time-varying rotation of the 
data constellation, which needs to be compensated for at the receiver 
before symbol recovery. This paper proposes a new approach for 
blind CFO estimation and symbol recovery. The received base-band 
signal is over-sampled, and its polyphase components are used to 
formulate a virtual Muhiplc-lnput Multiple-Output (MIMO) prob- 
lem. By applying blind MIMO system estimation techniques, the 
system response is estimated and used to subsequently transform 
the multiple CFOs estimation problem into many independent sin- 
gle CFO estimation problems. Furthermore, an initial estimate of 
the CFO is obtained from the phase of the MIMO system response. 
The Cramer-Rao Lower bound is also derived, and the large sample 
performance of the proposed estimator is compared to the bound. 

keywords-MuMi-user Systems, Distributed Antenna Systems, 
Carrier Frequency Offset, Blind MIMO System Identification 

1. INTRODUCTION 

In any communication system, the received signal is corrupted by 
carrier-frequency offsets (CFOs) due to the Doppler shift and/or lo- 
cal oscillators drift. The CFO causes a frequency shift and a time- 
varying rotation of the data symbols, which need to be compensated 
for at the receiver before symbol recovery. This can be achieved via 
pilot symbols. However, in the case of mobile systems and rich scat- 
tering enviroimients the effects of CFO become time varying and 
even small errors in the CFO estimate tend to cause large data re- 
covery errors. This necessitates transmission of pilots rather often, a 
process that lowers data throughput. In this paper we deal with CFO 
estimation without the need for pilot symbols. In single user sys- 
tems, or in multi-antenna systems in which the transmitters are phys- 
ically coimected to the same oscillator, there is only one CFO that 
needs to be estimated. This is typically done via a decision feedback 
Phase Lock Loop (PLL) at the receiver end. The PLL is a closed- 
loop feedback control system that uses knowledge of the transmitted 
constellation to adaptively track both the frequency and phase off- 
set between the equalized signal and the known signal constellation. 
However, depending on the constellation used during transmission, 
the PLL can have an M-fold symmetric ambiguity, resulting in a lim- 
ited CFO acquisition range, i.e., \Fk \ < l/(8Ts) for 4Q AM signals, 
where Ts is the symbol period. Moreover, the PLL typically requires 
a long convergence time. Alternatively, several methods have been 
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proposed [2], [4], [5], [7] [9] that blindly estimate the CFOs and 
recover the transmitted symbols using second-order cyclic statistics 
of the over-sampled received signal. Blind CFO estimation has also 
been studied in the context of orthogonal frequency-division multi- 
plexing (OFDM) systems, where the CFO destroys the orthogonality 
between the carriers (see [3] and the references therein). 

In a spatially distributed antenna system where data are transmit- 
ted simultaneously from multiple antennas, the received signal con- 
tains multiple CFOs, one for each transmit antenna. A PLL does not 
work in that case as there is no single frequency to lock to. The liter- 
ature on estimation of multiple CFOs is rather sparse. Existing liter- 
ature on this topic focuses on pilot based CFO estimation. In [6], the 
multiple CFOs were estimated by using pilots that were uncorrelated 
between the different users. In [13], multiple CFOs were estimated 
via Maximum Likelihood based on specially designed pilots. To ac- 
count for multiple offsets, [8] proposed that multiple nodes transmit 
the same copy of the data with an artificial delay at each node. The 
resulting system was modeled as a convolutive single-input/single- 
output (SISO) system with time-varying system response caused by 
the multiple CFOs. A minimum mean-square error (MMSE) deci- 
sion feedback equalizer was used to track and equalize the channel 
and to recover the input data. Training symbols were required in 
order to obtain a channel estimate, which was used to initialize the 
equalizer. 

Here we propose an approach for blind identification of multi- 
ple CFOs and subsequent symbol recovery. The received base-band 
signal is over-sampled, and its polyphase components are viewed 
as the outputs of a virtual MIMO system. The time-varying con- 
tribution of the CFOs, together with the transmitted symbols form 
the multiple-input/multiple-output (MIMO) inputs, while the time- 
invariant contribution of the CFOs along with fading channels com- 
prise the system response. By applying blind MIMO system estima- 
tion techniques, the system response is estimated and used to subse- 
quently transform the multiple CFOs estimation problem into many 
independent single CFO estimation problems. Furthermore, an ini- 
tial estimate of the CFO obtained from the phase of the MIMO sys- 
tem response can be used to initialize the PLL eliminating the sym- 
metrical ambiguity problem. The resulting method has full acquisi- 
tion range for normalized CFOs, i.e., \Fk\ < P/ {2Ts), where P is 
the over-sampling factor. To evaluate large sample performance we 
establish the Cramer-Rao Botmd (CRB) and compare the obtained 
mean square error to it. 

2. SYSTEM MODEL 

We consider a distributed antenna system, where K users trans- 
mit simultaneously to a base station. Narrow-band transmission is 
assumed here, where the channel between any user and the base 
station is frequency non-selective. In addition, quasi-static fading 
is assumed, i.e., the channel gains remain fixed during the packet 



length. The continuous-time base-band received signal y{i) can be 
expressed as 

K 

y{t) = auxk{t - Tfc)e^''"^'=* + w{t) (1) 
fc=i 

where Ofc represents the effect of channel fading between the fe— th 
user and the base station and also contains the corresponding phase 
offset; Tk is the delay associated with the path between the fc— th user 
and the base station; Fk is the frequency offset of the A;— th user and 
w{t) represents noise; Xkit) denotes the transmitted signal of user 
k: Xkit) = Sk{i)p(t — iTs) where Sk{i) is the i— th symbol of 
user fe; Ts is the symbol period; and p{t) is a pulse function with 
support [0, Ts]- 

Our objective is to obtain an estimate of s(i) = 

[si{i), SKii)]^ of the form 

s(i) = AP^s(i) (2) 

where P is a column permutation matrix and A is a constant diag- 
onal matrix. These are considered to be trivial ambiguities, and are 
typical in any blind inference problem. 

The received signal y{t) is sampled at rate 1/T = P/Ts, where 
the over-sampling factor P > if is an integer. In order to guar- 
antee that all the users' pulses overlap at the sampling times, the 
over-sampling period should satisfy: Ts/P > Tk,k — 1,...K, 
which means that the over-sampling factor P is upper bounded by 
Ts/min{Ti, tk}. I^et t = iTs + rriT, in = 1, P — 1 de- 
note the sampling times. The over-sampled signal can be expressed 
as 

ym{i) = viiTs + mT) 
fc=i 

K 

= Y. + + ^)'^ = 1' P - 1(3) 

where fk = FkTs/P, (\fk\ < 0.5) is the normalized frequency 
offset between the fc— th user and the base station, and the element 
of the virtual MIMO channel matrix A is given as 

am,k=ake^^'^"'^''pCpTs-Tk), m=l,...,P (4) 

Defining y(i) = [yi{i), ...,yp{i)]'^; A = {am,k}, a tall matrix 

of dimension P X K;s{i) = [si{i)e^^''^^'^ , ...,SK{i)e''^''^''^'^f; 

and w(i) = [w{i + -p), w{i + -p)]^, eq. (3) can be written in 
matrix form as 

y(i) = Ai(i) + w(i) (5) 

3. BLIND CHANNEL ESTIMATION AND 
COMPENSATION OF THE CEOS 

Let us make the following assumptions. 

• Al) For each m = 1...P, Wmi-) is a zero-mean Gaussian 
stationary random processes with variance crj, and is inde- 
pendent of the inputs. 

• A2) For each fc, Sfe(.) are a zero mean, independent identi- 
cally distributed (i.i.d.) stationary with nonzero kurtosis, i.e., 
7sj^ =C\im[sk{i),sl{i),Sk{i),sl{i)] 7^ 0. The s^'s are mu- 
tually independent, we can further assume that every user has 
unit transmission power, then Cs = I. 



• A3) The over-sampling factor P is no less than K. 

Under assumption (A2), it is easy to verify that tiie rotated input 
signals Sfc(.) are also zero mean, i.i.d, wide sense stationary with 
nonzero kurtosis. Also, the Sk(i)'s are mutually independent for 
different fe's. Assumption (A3) guarantees that the virtual MIMO 
channel matrix A in (5) has full rank with probability one. If the de- 
lays of users are randomly disttibuted in the interval [0, Tg/P), then 
each row of the channel matrix can be viewed as drawn randomly 
from a continuous distribution, thus the channel matrix has full rank 
with probability one. 

One can apply any blind source separation algorithm (e.g., [1]) 
to obtain 

A = APA (6) 

Subsequently, using a least-squares equalizer we can get an esti- 
mate of the de-coupled signals s{i), within permutation and diagonal 
scalar ambiguities as 

§(i) = (A^A)-^A'^y(i) = e'^''o^-^''\A\-'^V'^s{i?) (7) 

Denoting by 6k the fc— th diagonal element of ArgjA}, the fc— th 
separated input signal can be expressed as 

Ifc(i) = sfc(i)e^(-«'=+2''^'=*^) (8) 

Based on (8), any single CFO blind estimation method could be ap- 
plied to recover the input signal. Those methods can benefit by a 
CFO estimate provided by the channel matrix estimate as follows. 
The phase of the estimated channel matrix A equals 

/ 27r/i-F0i ... 2-KfK + ct>K \ 
* = Argk = P (9) 

V 2TvflP + ct>l ... 2-KfKP + <t>K j 

where </>fc = ^rg{afe}-|-6'fe, which accounts for botii the phase of afc 
and the estimated phase ambiguity in (8). The least squares estimate 
of fk is given by 

^'"2. P{Y.Up')-iT.Upy 

We can write fk = fk+ ffc where ek represents the estimation error. 

Noting that the de-coupled signals Sj{i) in (8) are shuffled in 
the same manner as the estimated CFOs in (10), we can use the es- 
timated CFOs to compensate for the effect of CFO in the decoupled 
signals (8) and get estimates of the input signals as 

s(i) = e^"'^''«^-'^>P^s(i) (11) 

Due to the residual error in the estimated CFOs, we can only 
compensate for a majority of the effect of CFO in (8) and obtain 

s^i) = s4i)e^(-^'=-^"^'=''^^ (12) 

Let us apply the PLL to the recovered signals Sj{i) in (12), to 
further mitigate the effect of residuary CFO ek . For 4QAM signals, 
as long as \P(-k\ < 1/8, it can be effectively removed by the PLL. 
Thus, the CFO estimator (10) can prevent the synmietric ambiguity 
of the PLL, and can also greatiy reduce the convergence time of PLL. 
From (9), we can see that the CFO estimator will achieve full acqui- 
sition range for the normalized CFO, i.e., \ fk\ < 1/2, which means 
we can deal with all continuous CFOs in the range Fk < P/{2Ts). 



4. CRAMER-RAO LOWER BOUND 

To evaluate the large sample performance of the proposed method, 
we establish the Cramer-Rao lower bound according to [12]. Via 
central limit theory arguments, the received signals y can be approx- 
imated as complex Gaussian signal, with zero mean, and covariance 
matrix given by 



AA^ + all 



(13) 



The covariance matrix is valid under assumption Al) and A2). The 
Gaussian assumption of the received signal is reasonable since the 
received signal is a linear mixture of i.i.d. signals. 
Let 



(14) 



where = [/i, /k]"^ is the vector of unknown CFOs, and 
= [ti, ...,rx]^ is the vector of random delays. The parameter 
we are interested in is the CFOs f , while p and cr^ are the nuisance 
parameters. 

Under the previous assumptions and the Gaussian approxima- 
tion, the Fisher Information Matrix (FIM) for the parameter vector 
a is given by [12] 



FIM, 



l,n = l,...,2K + 1 
(15) 

Since we are only interested in the CFG parameter f , following the 
derivation in [12], we can obtain that 



-!-CRB"^(f) = G^^G - G^"^A(A^A)"^A^G = G^n^G 

^ (16) 
where G and A are defined as 



lFIM = (^)-(C^«C-)(g|) 



A" 



[G A] 



where Cy = vsc{Cy) is a x 1 vector constructed from colimms 
of Cy, and G is of dimension x K, while A is of dimension 

X (K+l). 

To proceed, we just need to evaluate the derivatives of Cy with 
respect to a. First consider dcy/df^ , it holds that 



|^=vec&=vec([0.. 



• dfc • • • 0] A"-|-A[0 • 



•dfc • 



•Of) 



with dk 



(afe [1, P]'^), where © is the Hadamard ma- 
trix product. 

Similarly, we can get dcy/dp^: 

^ = vec(^) = vec([0 • • • e, ■ • • 0]A«+A[0 • • ■ ef • • • 0]^) 

withefe = [e2-4^MSzIt),...,e2-/^M2Wi)lT. 

Finally, we have that dcy/da'i, — vec(Cj^^). Now we have all 
the ingredients to evaluate CRB{f) from (16). 



delays, Tk, k = 1, .... K are uniformly distributed in the range of 
[0, Ts/P). The input signals used here are 4QAM signals. The es- 
timation results are averaged over 300 independent channels, and 20 
Monte-Carlo runs for each chaimel. 

The blind source separation algorithm used is the JADE 
method, which was downloaded from http://www.tsi.enst.fr/ car- 
doso/guidesepsou.html. 

We show the performance of both the pilots-based method and 
the proposed method at different data lengths and SNR set to 30dB. 
For the pilots method, each user transmitted a pilot signal of length 
32, and the pilots were random sequences uncorrelated between dif- 
ferent users. In Fig. 1 we show the Mean Squares Error (MSE) 
for the CFG estimator (10) for different values of the over-sampling 
factor P. To make the comparison fair for different over-sampling 
factor P, the MSE is calculated based on i ^^Li [(A - fk)Pf = 
1( X^^iK^fc ~ Fk)Ts\^ . We can see that by increasing P we can 
get more accurate estimates of the CFOs. In Fig. 2, we show the 
Bit Error Rate (BER) for different P's. For both the blind and the 
training based methods, the BER is calculated based on the recov- 
ered signals after the PLL. As expected, the BER performance also 
improves by increasing P. The proposed method appears to work 
well even for short data length. 

Next we show the performance of both methods at various noise 
levels. We set the packet length A*' to 1024. In Fig. 3, we show the 
MSE of the blind CFG estimator (10) as well as the training based 
method. We can see that by increasing P we can get more accurate 
estimates of the CFOs. In Fig. 4, we show the BER performance 
after PLL for both blind and training based methods. We can see 
that the proposed blind method has almost the same performance 
to the training based method for SNR lower than 20dB, while the 
training based method can achieve better BER performance for high 
SNR. 

The mean square error of the CFG estimator of (10) is plotted 
against the stochastic CRB derived in Section 4. In fig. 5, we plotted 
the MSE of the CFOs, as well as the CRB, as a function of the packet 
length T. We can see that the MSE curves are parallel to the CRB, 
and no error floor presented in the plot. Hence there is no bias in the 
estimates and the gap is only due to excess variance in the estimates. 
One possible reason for the existence of the excess variance is that 
we assume that we know the exact channel structure in the derivation 
of the CRB, i.e., the waveform used in transmission, which reduces 
the number of the unknown parameters. In the simulations, however, 
we did not assume any extra knowledge of the channel structure. 

We should note that the PLL is important for good symbol re- 
covery. For example, without the PLL, even if the residual error 
Ptk = Fk — Fk is only 0.001, the constellation will be rotated to a 
wrong position after 0.25/0.001 = 250 samples for 4QAM signals. 
To make sure that the PLL does not have the symmetrical ambiguity, 
we need to guarantee that I Pe/c I = \{Fk-Fk)Ts\ < l/8for4QAM 
transmission. Thus, on the average, the maximum tolerable MSE 
for the CFG is in the order of 10~^. From the simulations, we can 
see that the achieved compensation is sufficient for practical systems 
and commonly used modulation schemes. 



5. SIMULATION RESULTS 



6. CONCLUSION 



In this section, we verify the validity of the proposed method via 
simulations, under the following assiunptions. The channel coef- 
ficients ttk, k = 1, are zero-mean Gaussian random vari- 
ables. The waveform p{.) used here is hamming window. The con- 
tinuous CFOs are randomly picked in the range [— jf^i W')- The 



In this paper we have proposed a novel blind approach for identifica- 
tion of a distributed multiuser antenna system with multiple CFOs. 
By over-sampling the received base-band signal, we have converted 
the mulplte-input/single-output (MISG) problem into a MIMG one. 
Blind MIMG system estimation yields the system response, and 



MIMO input recovery yields the decoupled transmitted signals, each 
one containing a CFO. By exploring the structure of the MIMO sys- 
tems response we can obtain a coarse estimate of the CFOs, which 
can be combined with a decision feedback PLL to compensate for 
the CFOs in the decoupled transmitted signals. The proposed blind 
method has full acquisition range for normalized CFOs. We have 
provided a Cramer-Rao bound (CRB) for the proposed blind CFO 
estimators. The analytical results have been validated via simula- 
tions. 

7. REFERENCES 

[1] J.F. Cardoso and A. Souloumiac, "Blind beamforming for non- 
Gaussian signals," lEE Proceedings: Radar and Signal Processing, 
Vol. 140, No. 6, pp. 362-370 Dec. 1993. 

[2] P. Ciblat, P. Loubaton, E. Serpedin and G.B. Giannakis, "Performance 
analysis of blind carrier frequency offset estimators for noncircular 
transmissions through frequency-selective channels," IEEE Trans. 
Signal Process., Vol. 50, No. 1, pp. 130 - 140, Jan. 2002. 

[3] T. Roman, S. Visuri and V. Koivunen, "Blind frequency synchroniza- 
tion in OFDM via diagonality criterion" IEEE Trans. Signal Process., 
Vol. 54, No. 8, pp. 3125 - 3135, Aug. 2006. 

[4] M. Ghogho, A. Swami and T. Durrani, "On blind carrier recovery 
in time-selective fading channels," in Proc. 33rd Asilomar Conf. Sig- 
nals,. Systems, Computers, Vol. 1, pp. 243-247, Pacific Grove, CA, 
Oct. 1999, 

[5] F. Gini and G. B. Giannakis, "Frequency offset and symbol timing 
recovery in flat-fading channels: A cyclostationary approach," IEEE 
Trans. Comnmn., Vol. 46, No. 3, pp. 400 - 411, Mar. 1998. 

[6] F. Prihoda, E. Garbarine and A. P. Petropulu, "Resolving wireless 
collisions in random access networks," Proc. 2006 Asilomar Conf. 
Signals, Systems and Computers, Pacific Grove, CA, Oct. 2006. 

[7] K. E. Scott and E. B. Olasz, "Simultaneous clock phase and frequency 
offset estimation," IEEE Trans. Commun., Vol. 43, No. 7, pp. 2263 - 
2270, July 1995. 

[8] D. Veronesi and D. L. Goeckel, "Multiple frequency offset compen- 
sation in cooperative wireless systems" Proc. IEEE Global Telecom- 
munications Conf., San Francisco, CA, Nov. 2006. 

[9] Y. Wang, P. Ciblat, E. Serpedin and P. Loubaton, "Performance anal- 
ysis of a class of nondata-aided frequency offset and symbol timing 
estimators for flat-fading channels" IEEE Trans. Signal Process., Vol. 
50, No. 9, pp. 2295- 2305, Sept. 2002. 

[10] T. Kaiser and J.M. Mendel, "Covariance of finite-sample cumulants 
in array Ppocessing," Proc. IEEE Signal Processing Workshop on 
Higher-Order Statistics, pp.306-310, Banff, Alta., Canada, July 1997. 

[11] J. F. Cardoso, "Perturbation of Joint Diagonalizers" Tech. 
Rep. 94D023, Telecom Palis, Signal Department, 1994. 
http://citeseerist.psu.edu/cardoso95perturbation.html 

[12] P. Stoica, E. G. Larsson and A. B. Gershman, "The stochastic CRB 
for array processing: A textbook derivation," IEEE Signal Process. 
Letters, Vol. 8, pp. 148-150, May 2001. 

[13] O Besson and P Stoica, "On parameter estimation of MIMO flat- 
fading channels with frequency offsets" IEEE Trans. Signal Process., 
Vol. 51, No. 3, pp. 602- 613, Mar. 2003. 



MSE of CFOs at SNR=30dB 



10" 



10" 



o 
o 



10" 



P=2, bind 
-^t^ P=2, pilots 

- e - P=4, blind 

- + - P=4, pilots 




500 1000 1500 2000 2500 3000 3500 4000 4500 
packet length N 



Fig. 1. MSE of CFOs vs TV for K=2, with SNR=30dB, 4QAM 
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Fig. 2. BER vs for K=2, with SNR=30dB, 4QAM 
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Fig. 3. MSE of CFOs vs SNR for K=2, 4QAM 
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Fig. 4. BER vs SNR for K=2, 4QAM, T=1024 
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